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Abstract 

The dynamics of hydrogen bonds among water molecules themselves and 
with the polar head groups (PHG) at a micellar surface have been in- 
vestigated by long molecular dynamics simulations. The lifetime of the 
hydrogen bond between a PHG and a water molecule is found to be much 
longer than that between any two water molecules, and is likely to be a 
general feature of hydrophilic surfaces of organized assemblies. Analyses 
of individual water trajectories suggest that water molecules can remain 
bound to the micellar surface for more than a hundred picosecond. The activation 
energy for such a transition from the bound to a free state for the water 
molecules is estimated to be about 3.5 kcal/mole. 
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The study of hydrogen bond dynamics has proven to be a very useful tool |], ||, 
to understand the origin of many fascinating dynamical properties of water which 
are due to its extended hydrogen bond network ||. However, the same properties 
exhibit quite different behavior for water molecules at the surfaces of self-organized 
assemblies and biological macromolecules 0, ||, [7], |9|. Recent time domain spec- 
troscopic measurements have shown that their dynamics is considerably slower than 
their counterparts in bulk water, sometimes slower by more than two orders of mag- 
nitude^. While the reorientation of water molecules and solvation of ions or dipoles 
in bulk water proceeds with an average time constant of less than a picosecond (ps), 
the same at protein surfaces gets extended to hundreds of picoseconds [JTD], |7 , |q, |9 . 
Such slow dynamics has been observed in proteins [|7], [|, microemulsions \% |11 , 
micelles || [I2|] , lipid vesicles and bilayers [|TSj . 

The slow dynamics of water molecules at heterogeneous surfaces seems to be 
universal and could be a collective effect originating from the surface and/or from 
the nature of the hydrogen bond network at the surface. The work of Cheng and 
Rossky [|14| have already shown that water molecules at the surface of a protein 
can be structurally different from those in the bulk. A phenomenological theory 
developed recently proposes to explain the emergence of slow dynamics in terms of 
a dynamical equilibrium between bound and free water molecules at the surface [15| . 
Recent atomistic molecular dynamics simulations of an aqueous micellar assembly, 
the cesium pentadecaffurooctanoate (CsPFO), have shown that the orientational 
correlation function of the water molecules near the micelle show a dramatic slowing 
down in its long time decay |]I6|| . Note that micelles are not only of great interest 
themselves, but they also mimic surfaces of many biological molecules. The presence 
of a compact hydrophobic core and a hydrophilic surface resembles proteins rather 
closely. 

In this study, we explore the two most important microscopic aspects of water 
dynamics near the micellar surface - the hydrogen bond lifetime of the water-PHG 
hydrogen bond and also the distance dependence of the dynamics of the water-water 
hydrogen bond. The best way to study this dynamics is to investigate the hydrogen 
bond time correlation functions, S'hb(^) and Chb(^) (defined below), introduced 
originally by Rapaport [17j and used by Chandler et al to study pure water and 
more recently by Chandra || to explore effects of ions in water on the lifetime of 
the hydrogen bond. The time correlation function S^ 110 ^) for the hydrogen bond 
between the polar head group and the water molecules on the surface is defined by 



S£b (*) = (^wPHG(0)# w PHG(t)>/(/>wPHG> (1) 

where the population variable /i w phg(£) is unity when a particular tagged hydro- 
gen bond between a water molecule and a polar head group is hydrogen bonded 
at time t according to an adopted definition, and zero, otherwise. On the other 
hand, if w pHc(£) = 1 if the tagged polar head group- water hydrogen bond remains 
continuously hydrogen bonded during the time duration t and zero otherwise. Thus, 
S' H { g HG ()f:) describes the lifetime of a tagged pair. The angular brackets denote averag- 
ing over initial time values and over all water-PHG pairs. We employed a geometric 
definition of the wPHG hydrogen bond such that a water molecule was assumed to 
be hydrogen bonded to a surfactant if the distance between the oxygen of the water 
molecule and the carbon of the headgroup was less than 4. 35 A, and that the oxygen 
of the water molecule and the oxygen of the headgroup was less than 3.5A. These 
distance criteria were obtained from the first minimum in the respective pair corre- 
lation functions. For water-water hydrogen bonds, we have used the same criteria 
as used by Chandra 0, which were found to hold good for water- water hydrogen 
bonds near the interface also. 
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In order to understand the structural relaxation of a tagged hydrogen bond, one 
also defines the following time correlation function, C^ UG {t), 



(*) = (/>wPHG(OKpHG(t))/(/>wPHG> (2) 

Unlike S , gB HG (t), this correlation function is expected to exhibit a longtime tail fl], ||. 

The surfactant in this simulation is pentadecafluorooctanoate, with cesium be- 
ing the counterion, commonly referred to as CsPFO. The CsPFO-H^O system has 
been well studied experimentally ]TS| and is regarded as a typical binary to exhibit 



micellization. The amphiphiles form disc shaped micelles, stable over an exten- 
sive range of concentration and temperature. Details of our MD simulations have 



been discussed elsewhere [|16|, |19 |. The molecular dynamics simulation was carried 
out in the NVT ensemble for an aggregate of 62 CsPFO molecules in 10,562 water 
molecules. The potential for water molecules is the extended simple point charge 
(SPC/E) model | |20[| , the counterions carry a unit positive charge, which is com- 



pensated by a +0.4e charge on the carbon of the octanoate headgroup and a — 0.7e 
charge on each of the oxygens of the headgroup [21|. The equations of motion were 
integrated with the RESPA scheme |E| using the PINY-MD program |23] with an 
outer timestep of 4 fs. The analyses reported here were carried out from different 
sections of a subsequent 3.5 ns trajectory. The S(t) and C(t) functions were obtained 
with time resolutions of 12 fs and 1 ps respectively. The results reported here are 
at a temperature of 300K. 

Figure 1 shows the time dependence of the correlation function S'hb HG (^), with 
the inset showing the S^(t) (i.e., water-water) correlation function in bulk water, 
for comparison. Note the lengthening of the time scale in the H-bonding with the 
PHG of the micelle. S^g 110 ^) can be fitted to a sum of three exponentials, with 
the longest time constant being equal to 9 ps which is to be compared with 0.97 ps 
in pure water. The amplitudes and time constants are given in Table 1 From 
these, one can obtain an average time constant of around 6.8 ps, which is thirteen 
times larger than its value for the water-water hydrogen bond in pure water. It 
is difficult to pinpoint the reason for the slowdown of S , hb(^) at the surface. We 
found that the water molecules at the surface form bridge hydrogen bonds involving 
PHGs of the nearest neighbor surfactant molecules which could stabilize the wPHG 
hydrogen bond. 

In figure 2, we show the hydrogen bond time correlation function Cjjb HG (^)- I n 
the inset we show C^(t) for bulk water. Note again the lengthening of the long 
time decay. The values of the time constants and the respective amplitudes are given 
in Table 1. In this case, the long time constant stretches to more than one hundred 
picoseconds. The reason for this unusual long decay time can be traced back to 
those trajectories which leave the mi cellar surface to go to the bulk but return after 
a long time to get bonded to the same PHG at the surface. This is in line with the 
explanation of the long time decay of the Chb(^) function in bulk water M, p5[]. In the 
present problem, the existence of such trajectories indicate presence ofcorrelations 
in the surface region. The S(t) function is thus a more accurate representation of the 
lifetime dynamics of the hydrogen bond than the C(t) function, although the latter 
contains rich information on the correlated pair diffusion of the water molecules. 

Analysis of the trajectories of the individual water molecules at the micellar 
surface reveals an amazing richness of events. In figure 3, we show representative 
trajectories over 800 ps, of four arbitrarily chosen water molecules that were hy- 
drogen bonded to the surface at time zero. The plot shows the shortest distance, 
-D\y-PHG) of each of the water molecules to the micellar surface as a function of 
time |26| . For example, the water molecule shown in the top left figure, rattles near 



the surface for about 150 ps, then diffuses into the bulk region, stays there for around 
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500 ps, and then revisits the micelle surface, presumably to form a new hydrogen 
bond with a PHG. It should also be noted that the micelle itself is quite fluxional 
over these time scales. Note also that while the Shb HG (*) and C^ KG (t) functions 
study the dynamics of a particular wPHG hydrogen bond (with average lifetimes of 
6.9 and 43.7ps respectively), the apparently longer lifetime of bound water molecules 
shown in the trajectories in figure 3 arises from the fact that a water molecule can 
remain bound to the micellar surface even after its hydrogen bond with a particular 
PHG is broken. 

Such dynamical processes involving water molecules affect the way they form 
hydrogen bond between themselves. We have explored this aspect by studying the 
water-water hydrogen bond time correlation functions, for water molecules that 
belong to different interfacial layers. In figure 4 we show the lifetime correlation 
functions, Cgg 7 (£) and S'hb W between two water molecules that both exist at several 
regions away from the surface p6| . This function shows a sharp slow down as the 



surface is approached closely. The time constant obtained from the Shb function for 
the water molecules within 6A from the micellar surface is around 25-30% larger 
than that for water molecules in the bulk region. The time constant for the same 
water molecules obtained from the Chb function is about 45% slower than that for 
bulk water. Again, the S'hb function decays much faster than the Chb one. 

The simulation results can be combined with a theoretical model to obtain useful 
information on the micelle-water interaction. The model assumes a dynamical equi- 
librium between the "bound" and "free" states of water molecules at the micellar 



surface |15|]. The dynamical variable h(t) represents the instantaneous population of 
the bound state and the correlation function, C^ RG (t) gives, under the regression 
hypothesis (or linear response theory) |27H , the decay of an initial (t=0) bound state 
and its formation at a later time. Therefere, we can set, in the long time, Cjjb (t) 

~ e t bf , where t^f is the time constant for the bound to free transition. From 
Table 1, we find the average time constant at 300K to be 43.7ps. We can obtain an 
activation energy E\, for this transition using transition state theory to express the 
rate constant (— ), as, 

\ TBF ' ' ' 

1 k B T -1a_ 
t B f h 

where, k-Q is the Boltzmann constant, T is temperature, and h is the Planck con- 
stant. E\ is then found to be 3.34 kcal/mole. A similar analysis of C H / B HG (t) at 
350K (data not shown here) yields a time constant of 18.9ps. Using an Arrhenius 
dependence of the time constant with temperature, we find that the activation en- 
ergy is 3.50 kcal/mole, in close agreement with the result of the transition state 
theory presented above, showing the robustness of the above analysis. 

To conclude, we have carried out large scale atomistic molecular dynamics simu- 
lations of an aqueous micellar system, with emphasis on understanding the hydrogen 
bond breaking dynamics at the micellar surface. We find the rather surprising result 
that the hydrogen bond between the micellar polar head group and a water molecule 
has a much longer lifetime - almost 13 times larger than that in the bulk between 
two tagged water molecules. This longer lifetime could originate from the bridge 
bonds that the surface water molecules form with the PHGs and also from the cou- 
pling to the micelle which acts as a bath for HB excitations. Another interesting 
result is that the lifetime of hydrogen bonds between two tagged water molecules 
increase only by about 25-30% as one approaches the micellar surface. Although 
this slowdown is not as dramatic as the slowdown in the orientational relaxation 
of the water molecules with large residence times near the interface, it is rather 
sharp. Analysis of the individual trajectories reveal that the quasi-bound water at 
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the micellar surface (immobilized by double or more hydrogen bond bridges) may be 
partly responsible for the slow dynamics. The activation energy for the transition 
from a bound to a free state for the water molecules on the micellar surface, esti- 
mated from both the transition state theory and from the temperature dependence 
of the hydrogen bond lifetime, is about 3.5 kcal/mole. The agreement between the 
rather simple model and the temperature dependent simulation data is remarkable. 
We believe that the results discussed here are generic to organized assemblies and 
biological macromolecules that possess a hydrophilic surface or hydrophilic pockets 
that are accessible to water. 
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Table 1: Parameters of multi-exponential fits to the water-PHG hydrogen bond time 
correlation functions shown in Figures 1 and 2. 



Function 


Time constant [ps] 


Amplitude (%) 




0.3 


5.9 


Shb(*) 


3.6 


31.5 




9.1 


62.6 




3.4 


16 


Chb(*) 


29.0 


63 




118.5 


21 
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Figure 1: Sgg HG (t) function for the hydrogen bond between the polar head group 
and water molecules. Inset shows the same between pairs of water molecules in pure 
water. 
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Figure 2: Cgg HG (t) function for the hydrogen bond between the polar head group 
and water molecules. Inset shows the same between pairs of water molecules in the 
region far away from the micelle [J2^ J. 
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Figure 3: Trajectory of four water molecules for a time period of 800 ps. The 
shortest distance of the water molecules to the micellar surface, -Dw-phg; is plotted 
against time. 
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Figure 4: The location dependence of the time correlation function, Cg^t), for the 
hydrogen bond between pairs of water molecules located at different regions from the 
mi cellar surface [p6| . Solid curve: Within 6A ; Long-Dashed curve: Between 6A and 
9A ; Dashed curve: Beyond 25A . Inset shows the similar location dependence of 
the Sgg(t) function for the hydrogen bond between pairs of water molecules. 
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